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Spontaneous global symmetry breaking of 0(3) scalar field gives rise to point-like topological 
defects, global monopoles. By taking into account self-gravity, the qualitative feature of the globa l 
monopole solutions depends on the vacuum expectation value v of the scalar field. When v < 
th ere a re global mono pole solutions which have a deficit solid angle defined at infinity. When 
y/l/87r < v < -\/3/87r, there are global monopole solutions with the cosmological horizon, which 
we call the supermassive global monopole. When v > ^3/8n, there is no nontrivial solution. It 
was shown that all of these solutions are stable against the spherical perturbations. In addition 
to the global monopole solutions, the de Sitter s olutio ns exist for any value of v. They are stable 
against the spherical perturbations when v < ^3/8%, while unstable for v > ^3/8n. We study 
polar perturbations of these solutions and find that all self-gravitating global monopoles are stable 
even against polar perturbations, independently of the existence of the cosmological horizon, while 
■ the de Sitter solutions are always unstable. 

O \ PACS numbers: 11.27.+d, 14.80.Hv, 04.40.-b, 98.80.Cq 
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I. INTRODUCTION 



Phase transitions in the early universe are caused by the symmetry breaking leading to a manifold of degenerate 
vacua with nontrivial topology and giving rise to topological defects. The defects are classified by the topology of 
the vacua such into domain walls, cosmic strings and monopoles. If the gauge field is involved in the spontaneous 
symmetry breaking, the defects are gauged. On the other hand, when the symmetry is global, the emerging defects 
■ are called global defects. 

^vq I In this paper, we shed light on global monopoles. Although energy of the gauge monopoles is finite, the global 
monopoles have divergent energy because of the long tail of the field. This divergence has to be removed by cutting 
""^ ' off at a certain distance. This procedure is not necessarily artificial, because another defect which may exist near the 

^jH, original one cancels the divergence. This secondary defect is not only the monopole, but also can be a domain wall 

^ or a cosmic string. 

Global monopoles have been an interesting subject in cosmology. They were thought as of the seeds of structure 
formation or inflation. By taking into account the self-gravity of the global monopoles, it can give rise to a deficit 
solid angle Q, which would affect cosmological data. It, however, may be counteracted when the universe has a 
cosmological constant. 

Vilenkin and Linde independently pointed out that topological defects can cause inflation |^,|| . When the vacuum 
expectation value (VEV) is larger than a certain critical value, the scalar field stays on the top of the potential after 
global symmetry breaking. In this case, the space would expand exponentially with time. This inflation model is 
called topological inflation. Topological inflation does not suffer from the initial value problem, unlike the new or 
hybrid inflation. Sakai et al. 0] found the critical vacuum expectation value (VEV) is 0.33-Mp numerically. 

Recently new types of the self-gravitating global monopole solutions were discovered numerically || . As the VEV 
of the 0(3) scalar field v increases, the deficit solid angle also gets large and it becomes 47r when v — v cr i := y/l/8n. 
Beyond this critical value there is no ordinary monopole solution, but there appears a new type of solution in the 
parameter range v cr t < v < v max := y/3/8n. This has a cosmological horizon at r = IL. and gives a natural cutoff 
scale. The appearance of the new solution is similar to the supermassive string solution p] , which has a deficit angle 
larger than 2n. In this sense, we will call this solution the supermassive global monopole in contrast to the ordinary 
global monopole. 

From gravitational theoretical interest, the black hole counterpart of the global monopole is also discussed. Maison 
investigated the scalar hair of the black holes in the same system and showed their existence domains as a function of 
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the radius of the event horizon M. Nucamendi and Sudarsky discussed the definition of the Arnowitt-Deser-Misner 
(ADM) mass in spacetime with deficit solid angle and found that it becomes negative for small black holes § . 

One of the most important issues of these kinds of isolated objects is stability. In the ordinary monopole case without 
gravity, if Derrick's no-go theorem M could be applied, they would be unstable towards radial rescaling of the field 
configuration. This is, however, not the case due to the diverging energy of the solutions. It was demonstrated that 
the ordinary monopole solutions are stable against spherical perturbations. As for the non-spherical perturbations, 
there was some debate. Goldhaber |[o| investigated the polar perturbations and found that the energy functional is 
the same form as that of the sine-Gordon equation under some conditions. Hence there would be a zero mode leading 
a north-pointing teardrop shape for the monopole mass density and the knot would be untied. Rhie and Bennett 
fj"lf , however, pointed out that such instability is just the artificial fixing of the monopole core. If the monopole is 
free to move, which is the natural situation of the isolated system, the monopole cannot become a teardrop shape 
and is stable. This was actually confirmed by numerical simulations ]l^Jl^| . Bennett and Rhie (lj] also pointed out 
by numerical calculation with 2D code that a monopole and antimonopolc pair would collapse to a string through 
the unwinding process when these cores are artificiality fixed. If these cores of the pair are free to move, such an 
unwinding process does not occur. 

It is expected that the self-gravitating global monopoles with v < v cr i (i.e., without cosmological horizon) have the 
same stability properties as a non-gravitating one. Stability may change, however, for the supermassive monopole 
solutions (v > v cr i), because the domain of communication, that is, the boundary condition, is different due to the 
cosmological horizon. Maison and Liebling investigated the spherical perturbations of the supermassive monopole 
solution [ jl5| . They make use of de Sitter solutions, which are trivial solutions such that the scalar field stays at the 
top of the potential barrier and exists for any value of the VEV of the scalar field. By their analysis the stability 
change of the de Sitter solutions occurs at v = v max , beyond which the solutions are stable, while unstable below 
that. And the supermassive monopole solutions emerge just at this value if the VEV decreases from a larger value. 
This kind of behavior can be seen in a variety of systems in nature and explained by using catastrophe theory. The 
supermassive monopole solutions inherit the stability from the de Sitter solution with v > v max . As a result, they are 
stable against the spherical perturbations even if they have a cosmological horizon. 

Then, are the supermassive monopole solutions really stable? We have to examine this question carefully. This is 
because the polar perturbation pushes the scalar field configuration of the de Sitter solution to one direction in the 
internal space from the top of the potential barrier. In the anti-de Sitter background, such a solution can be stable 
even in the tachyonic situation if the effective mass of the scalar field satisfies the Breitenlohner-Freedman bound jl6) . 
In the de Sitter case, however, it is easy to imagine the scalar field rolls down to its VEV, i.e., the solution is unstable. 
Hence, the supermassive monopole solutions may be unstable against the polar perturbations. To settle this issue is 
the main purpose of this paper. 

This paper is organized as follows. In Sec. |l|we review the static solutions and their spherical stability. In Sec. IH 



we show the instability of the de Sitter solution against the polar perturbations. In Sec. IV, wc formulate the polar 
perturbations of the 0(3) scalar field. In Sec. [V], we show the stability of the self-gravitating global monopoles. 
Throughout this paper, we use the units h = c = G = 1. 



II. STATIC SOLUTIONS 



In this section, we briefly review the self-gravitating global monopole solutions Jl|,fl5[]. The theory of a scalar field 
with spontaneously broken internal 0(3) symmetry, minimally coupled to gravity, is described by the action 



S = Id x\J—g 



R 
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A 



(1) 



where R is the Ricci scalar of the spacetime and $ a (a = 1, 2, 3) is the triplet scalar field. A and v are the self-coupling 
constant and the VEV of the scalar field, respectively. The energy momentum tensor is 



^d p $ a d p <5> a + ^($ 2 -v 2 ) 2 



For the static solution with unit winding number, we adopt the so-called hedgehog ansatz 

= h{r)— , 



(2) 



(3) 



where x a are the Cartesian coordinates. 
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We shall consider the static spherically symmetric spacetime and adopt a Schwarzschild type metric 



ds 2 = -f{r)e- 2S(r) dt 2 + Jj^dr 2 + r 2 (dd 2 + sin 2 6dcj) 2 ), 



where 



m = i - 



2m(r) 



Under these Ansatze, we get these field equations 
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m' = Airr 
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X(h 2 - v 2 )h, 



(8) 



where a prime denotes a derivative with respect to the radial coordinate. 

These equations are integrated with suitable boundary conditions. At the center the spacetime should be regular. 
By expanding these equations, we find h'(0) can be regarded as a free parameter, which is determined by the other 
boundary condition at r — » oo (for the ordinary global monopole case) or r = R c (for the supermassive global 
monopole case). For the ordinary global monopole solution the spacetime approaches asymptotically flat spacetime 
(which implies that the curvature vanishes) with deficit solid angle a 
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On the other hand, we impose the existence of the regular cosmological horizon at r — R c for the supermassive global 
monopole. 

There is a trivial de Sitter solution h(r) = 0, f(r) = 1 — r 2 /R 2 , and S(r) = 0. This solution has a cosmological 
horizon at r — R c — ^/3/A e //, where the effective cosmological constant is A e // := 2tt\v 4 . These solutions exist for 
any value of v. 

The ordinary global monopole solutions exist for < v < v cr i. The configuration of the scalar field is shown in 
Fig. |l| (v = 0.15). We set A = 0.1 (which is adopted in Ref. ||) without loss of generality throughout this paper since 
A can be scaled out by introducing new variables f := A 1 / 2 wr, fh := \ x l 2 vm and i> a := Q a /v. The deficit solid angle 
a becomes large as a = 47r(87ru 2 ), and a — Air for v — v cr i, which implies the disappearance of the asymptotic region. 

Beyond the critical value the supermassive global monopole solutions appear for v cr i < v < v max 0- This has 
a cosmological horizon. If y / 2/87r < v < v ma x, the scalar field shows the oscillating behavior approaching its VEV 
beyond the cosmological horizon (See Fig. |l|). Asymptotically it becomes the de Sitter spacetime pl|. At v = v max 
the solution coincides continuously (at least in the domain of the communications) with the de Sitter solution. This 
can be seen in the behavior of the parameter h'(0) as shown in Fig. ||. 

Maison and Liebling investigated the stability of the de Sitter solutions against spherical (both in spacetime and 
in internal space) perturbations and found that stability changes at v = v max . By this result they expected that the 
stable property of the de Sitter solution with v > v max is transferred to the supermassive global monopole solutions. 
This kind of study was performed by using catastrophe theory and applied for the black hole spacetimes |p7| , |l8| . 
However, are the supermassive global monopoles really stable even for the non-spherical perturbations? The 0(3) 
field which constructs the de Sitter solution is always on the top of the Mexican hat potential, which seems unstable. 
It is easily imagined that the scalar field rolls down to their VEV by just pushing it in one direction in the internal 
space. We investigate this problem in the next section. 
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FIG. 1. The field configurations (above) and the metric function f(r) (below) of the global monopole solutions for v = 0.15 
(dashed line), 0.25 (solid line), 0.30 (dot-dashed line) and 0.34 (dotted line). The solutions except for v = 0.15 have a 
cosmological horizon. The solutions for v > y/2/8n w 0.28 show oscillating behavior. 




FIG. 2. The diagram of v versus h! at r = 0. The solid and the dashed line represent the monopole and the de Sitter solutions, 
respectively. The stability here means against the spherical perturbations. This figure says that the monopole solutions and 
de Sitter solution with v < v max are stable against the spherical perturbations. v max is the bifurcation point which connects 
the stable de Sitter, the global monopole, and the unstable de Sitter solutions. 



III. INSTABILITY OF DE SITTER SOLUTIONS 



Now, we investigate time dependent perturbations of the de Sitter solutions which are non-spherical in the internal 
space. We perturb one component of the 0(3) scalar field 



^(t,r) = $ 2 (t,r) ee 0, $ 3 (i,r) 



Acrt 



C(r). 



(12) 



Note that the metric functions are not affected by this perturbation since the back reaction is second order, and we 
can consider the de Sitter spacetime as a background. If the perturbation equation allows a solution with imaginary 
a, the mode function evolves exponentially with time. That indicates instability of the de Sitter solution. 
The perturbation equation is written as, 



1 



[r 2 /C']' + A« 2 C = -/-VC. 



(13) 



Adopting the tortoise coordinate r*, 



dr 
dr* 



(14) 



and a new variable £ = r£, we can rewrite it in the Schrodinger equation, 



(15) 
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U(r) is the potential function of the linear equation 



U(r) = -f\v 2 (l + ^ttv 2 ^ . (16) 

The potential function never becomes positive inside the cosmological horizon. This is the exactly the same form as 
the Eqs. (37)-(39) in Ref. fl|] in the de Sitter case. In that paper it was shown that these equations always have at 
lease one negative eigenmode for any value of A and v. Hence, it is concluded that all the de Sitter solutions which 
we consider here are unstable. 



IV. POLAR PERTURBATION OF THE SCALAR FIELD 



Since the de Sitter solutions have polar instabilities, the self-gravitating global monopoles may inherit them. Hence 
we study the polar perturbation of the global monopoles in general relativity carefully in this section and the following 
section. 

The polar deformation of the 0(3) scalar field was proposed by Goldhaber JhJ. He introduced a new coordinate 
y := lntan(0/2) to discuss the invariance of the energy under the deformation. Achucarro and Urrestilla improved 
his notation Q, 



4- 1 = H(t, r, 9) s'md(t, r, 9) cos </>, 
$ 2 = H(t, r, 9) sin S(t, r, 9) sin <j>, 
$ 3 = H(t, r,9) cos 0(t,r,9), 



(17) 



and 



tan(0/2) 



_ p y+£(t,r,9) 



(18) 



is the polar component of $ a and H (t, r, 9) = h(r) + 6h(t, r, 9). When £ = 0, i.e., 9 — 9 and Sh — 5h(t, r) the global 



monopole is spherically symmetric, while it become a 'string' when £ 
is expressed with new coordinate y, 



oo. The energy of the static global monopole 



E = drdyd(j)(pi + p 2 ), 



(19) 



where 



Pi 



IP 
2 

1 



cosh 2 y 



P2 = ~ 2 



dH 

dy 



2 cosh y 



dH 
dr 



,,2\2 



(20) 



(21) 



In the far region from the monopole core, the terms dH/dr and dH/dy can be disregarded. Goldhaber pointed out 
that the first two terms in p\ are in the same form as the energy of a sine-Gordon soliton, which implies that the 
energy is invariant under translation of coordinate y (or £). Thus, he concluded that the global monopoles have 
instability if there is a deviation in which H = v and 89 /dr — are held. Rhie and Bennett proved, however, that 
such a collapse does not occur, if the monopole core is free to move JyJ. 
We assume £ and Sh are small. Therefore we take 



9 = 9 + £ sin 9. 



(22) 



Now we have two perturbative functions £ and Sh for the 0(3) scalar field. The simplest polar perturbation assumes 
that Sh = and £ is independent of 9. However, this is inconsistent. In a non-relativistic case, such perturbation 
must vanish through a perturbation equation. If £ has 9 dependence, we will get 



<9 e £ + 2£cot6< = 0. 



(23) 
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FIG. 3. The eigenmodes of rj and h£/v for a = 0.1 (solid line) and 0.001 (dashed line) when v = 0.25. We can see the 
oscillating behavior. These are continuous modes. 

The solution of this equation is £ oc (sin#)~ 2 , which is unphysical because £ violates the perturbative approximation 
around the axis. 

The same feature appears in the case that only Sh is considered. Consequently, we can conclude that perturbation 
only with either £ or Sh does not occur when the self-gravity is set to zero. We will find that this is true also in the 
self-gravitating case. 

Now, we assume £ is independent of 9. Since we will examine the I — 1 Legendre type perturbation below, we take 
Sh = rj(r) cosde t<Tt . The perturbed scalar field is 

(& 1 = (h + r\ cos 9e l(7t ) sin 9 cos 4> + h£_e l<7t sin 9 cos 9 cos <fi, 
<J> 2 = (h + 7]cos9e l(Tt ) sin 9 sin (f> + h^e l(Tt sin 9 cos 9 sin (f>, 
$ 3 = (h + rj cos 9e lat ) cos 6 - hie lat sin 2 9. 

(24) 

Let us consider the combination 

p(r) = h£ + Ti, q{r) = h£. (25) 
in Eq. (p4|). Putting p(r) = 0, we recover the polar perturbation of the de Sitter solution 

$ 1 = $ 2 = 0, $ 3 = -q{r)e l,Tt . (26) 



which is discussed in Sec. [II. 



V. POLAR PERTURBATION OF THE GLOBAL MONOPOLE 

Achucarro and Urrestilla studied stability of the global monopoles against polar perturbation with all orders []1J 
and found that the global monopoles are stable. Here we extend their analysis to the self-gravitating cases. 

For the global monopole solutions, the matter field is non-zero. So the first order perturbations of the matter field 
couple to 0-th order to create the metric perturbations, which can be described as 

ds 2 = -fe- 25 e 5v dt 2 + je s ^dr 2 + r 2 e s ^d9 2 + r 2 sin 2 9e^dcf> 2 . (27) 

The valuables in the metric can be separated , 

6v = ^2Ni(t,T)Pi(co89), (28) 
Sna=^2Lt(t,r)Pi(coB0), (29) 

5^3 = J2 t T * (*' r ) P * ( cos ) + Sl (*' r ) P ^.« ( cos *)] ' ( 3 °) 

6fl> = t T ' (*» r ) p i ( cos e ) + S i (*« r ) p ifi ( cos °) cot °i ( 31 ) 



G 



where Pi is a Legendre polynomial. Now we consider the simplest case I — 1 and drop the suffix /. For the higher 
order Legendre polynomials, the eigenvalue a 2 is expected to be larger than that for ( = 1. Thus, the perturbed 
metric can be written in the form, 

ds 2 = - [/ + B{r)e lat cos<9] e- 2S dt 2 + j [l + L(r)e i<Tt cos(9] dr 2 + r 2 [l + T{r)e i<Jt cosd] (d6 2 + sin 2 9d<f) 2 ). (32) 

Here, we introduced a new variable B(r) := N(r)f(r) for convenience, re-defincd T := T — S and assumed harmonic 
time dependence. 

The energy momentum tensor is, 



T t 


(0) 

= t\ + 


fh' 2 2h h 2 


- Xh(h 2 


-v 2 )r] cos0e i<rt 


T/ 


<°) r 

= Tr + 




\h(h 2 - 


v 2 )r) 


cosfle" 7 *, 




= in, ftirj cos de la \ 










h 2 









(33) 
(34) 

(35) 
(36) 



(0) (0) 

where T \ and T r r are non-perturbed components. We displayed only the first order perturbation. Thus, we can get 
the perturbation equations 



L + T= 16tt/i 2 £, 



f(-rT' + L- rTS' - T) + -f'rT = 8nrfh'r], 



(37) 
(38) 



fr 2 T" - {frL)' + -f'T'r 2 + 3/rT' - L = 4tt [-2fr 2 h'r]' + fr 2 h' 2 L + 2h 2 T - 4h(h£ + rf) -2r 2 Xh(h 2 - v 2 ) V ] , (39) 



fr 2 T" - frL' -L + 2frT' - 2frS'L + r 2 fS'T' 
from the Einstein equation, and 



r 2 a 2 e 2S 



T-B'r+ — + 



B Brf 



f f 



Snr 2 fti (h'L + 2rj') . 



a 2 he 2S 
f 



s + m)" + i-f+f'-s'f\(hsy- 



277 + 2/i£ hL 1 hB 



(40) 



(41) 



a 2 e 2S 



V + fv"+[-f + f'-S'f)v' + fh'T'- 



2h 



+ Xh(h 2 -v 2 ) 



L — 



4/i£ + 4rj- 2hT 



fh'L' 1. , h'f'B w i2 2 , 



(42) 



from the equation of the scalar field. 

The boundary conditions at r = are obtained by imposing regularity. By Taylor expansion around r = 0, we find 



£ = &r + O (r 3 ) , 
hi 



(86 - BOr 2 + O (r 4 ) , 



S = Sir - 1 (m 3 + Y nh i) B ^ + ° ( r5 ) ' 

T = InhjBir 3 + O (r 5 ) , 
5 

L = (l6n^ - \vB^j hjr 3 + O (r 5 ) , 



(43) 
(44) 

(45) 

(46) 

(47) 
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FIG. 4. The diagram of the shooting parameter Bi versus a 2 when v — 0.15. The background global monopole solution does 
not have a cosmological horizon. We can see a 2 converges to positive constant for large B±. The horizontal line is minimum 
eigenvalue a m i„ := 2Xv 2 in the non-self-gravitating case. The asymptotic value of a 2 for the self-gravitating case is smaller 
than that of the non-self-gravitating case. 



where F n represents the n-th order derivative coefficient of F at r — 0. We assumed all values at r = are 0, for the 
elimination of trivial translations mode in which rj = — h' cos9 and h£ = h/r Jl3| . This assumption is not, however, 
an artificial fixing of the core, but only coordinate transformation. The values of hi and 7713 are given by the static 
solutions in Sec. |J. £1 an d B\ cannot be determined by this regularity condition. Among these, £1 is arbitrary 
because of the freedom of the constant multiplication in the linear theory. Hence there are two parameters B\ and 
tr in this system, which should be adjusted to obtain the regular normalizable eigenmodes. If there is no solution, 
we cannot find an appropriate set of B\ and a 2 . 

By these boundary conditions we solve Eqs. (|37])-([l2]) numerically. Fig. [5] is the typical solution of the supermassive 
case. We can find the oscillating behavior for the positive eigenvalue a 2 . These are the continuum modes. 

If we neglect self-gravity, the perturbation equations become quite simple. In this case it is easy to observe that the 
potential function is nonzero constant asymptotically to infinity. Hence there is a minimum value of the eigenvalue 
a min = 2^« 2 > for the continuous modes. When self-gravity is taken into account, the situation does not change 
seriously for the ordinary global monopole solution. The existence of the minimum eigenvalue is seen in Fig. [|. 
Although o m in decreases continuously for large B\ , it converges to a positive constant which is smaller than that of 
the non-self-gravitating case. 

For the supermassive global monopole, however, the form of the potential function is qualitatively different due to 
the existence of the cosmological horizon. It vanishes at the horizon without bottom up, and hence the continuous 
modes exist for the infinitesimally small eigenvalue as seen in Fig. [| 

Figure ^ shows the "zero modes" (a 2 = 0), which diverge and are non-normalizable, for several values of v. We 
cannot find a real zero- mode solution for any value of B\ while there exists a solution for non-zero positive a 2 . When 
a negative eigenmode exists in this kind of linear perturbation analysis, the perturbed functions have at least one zero 
point for a 2 = in general. In our case, "zero modes" are positive everywhere. This indicates that the self-gravitating 
supermassive global monopoles are stable against the polar perturbations as well as the spherical perturbations. As v 
approaches its maximum value v max , the "zero modes" become small. It is expected that it has zero when v = v max . 



This behavior is consistent with the instability of the de Sitter solution discussed in Sec. Ill 



VI. DISCUSSION 



We investigated the polar stability of the self-gravitating supermassive global monopoles and de Sitter solutions in 
the Einstein-0(3) scalar system by the linear perturbation method Although the de Sitter solutions always have 
at least one unstable mode for any value of the VEV of the scalar field, the supermassive global monopole solutions 
do not. This implies that the supermassive global monopoles are stable against polar perturbations. We also find that 
the minimum eigenvalue is infinitesimal for the supermassive global monopoles while it becomes non-zero finite for the 
ordinary self-gravitating global monopoles and non-gravitating counterparts. This is due to the different boundary 
conditions at the cosmological horizon. 

Our analysis can be extended to the black hole solution inside of the monopole, i.e., the black hole solution with 
0(3) scalar hair. Such solution was discovered a decade ago and its supermassive counterpart was recently 

discovered by Maison and Liebling 111] . It was reported that these solutions are stable against spherical perturbations. 
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FIG. 5. Diagram of the shooting parameter B\ versus a 2 when v = 0.25, which corresponds to the case of the supermassive 
global monopoles. We can see that a 2 is continuously decreasing to 0, but there is no zero mode as seen in Fig. o. 




FIG. 6. "Zero modes" of the supermassive global monopoles. We can see that the "zero modes" are positive and have no zero 
point. This indicates the stability of the supermassive monopoles. At v = ^3/8% w 0.345, the supermassive global monopole 
coincides with the de Sitter solution, which is unstable. This figure shows that the "zero modes" would have zero point as 

V ► Umax • 



This result is notable because this scalar hair can be a physical hair. Hence we should check whether or not this hair 
is stable also against polar perturbations. It should be noted however, that even if this hair is absolutely stable, it 
does not mean the violation of the the black hole no-hair conjecture, because such conjecture is implicitly assumed in 
the asymptotic flatness. Some attempts to extend the conjecture to asymptotically non-flat spacetime were discussed 
in Ref. 
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